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I. INTRODUCTION 

Since its theoretical prediction in 1962 Q, the Josephson effect has played a major role in the physics and technology 
of superconductors The physics of the Josephson effect becomes manifest when a weak link is created (e.g. by 
tunneling or through a point contact) between two systems which have undergone some type of gauge symmetry 
breaking. In fact, the Josephson effect has also been observed in superfluids ||, and, because of the profound 
analogies, its observation in the recently achieved Bose condensed atomic gases is generally expected. A system 
whose gauge symmetry has been spontaneously broken can be described by an order parameter that behaves in many 
respects like a macroscopic wave function ^(r). In the simplest cases, the order parameter reduces to a complex 
scalar, ^> = ^fpe 1 ^ , where p(r) is the "superfluid density" and (f>(r) is the phase. For Bose-Einstein condensates of 
dilute alkali gases, ^(r) is the wave function of the macroscopically occupied one-atom state. 

In his epoch making paper JjJ, Josephson predicted that, between two weakly connected superconductors of phases 
ipi and (f2, a non-dissipative particle (Cooper pair) current flows between them whose value is 

I (if) = I c sin <p, (1) 

where I c is the critical current and tp — ip2 — tpi is the relative phase [<j>(r) is assumed to be approximately uniform 
within each superconductor]. He also predicted that, in the presence of a nonzero chemical potential difference 
fj, = 112 — fJ-1, the relative phase rotates as 

<P = (2) 

The Josephson relations can be obtained from very general considerations and apply to any pair of weakly linked 
systems that can be described by macroscopic wave functions Q] . They can be obtained as the equations of motion 
of the "pendulum Hamiltonian" 

H(<p ) N) = Ej(l-cos<p) + ±E c N 2 > (3) 

where Ej = hl c is the Josephson coupling energy, N — (N2 — N\)/2 is the number of transferred particles, and 
Ec = dfi/dN is the capacitive energy due to interactions. In the absence of external constraints, /1 = EqN. In a 
superconducting link, the critical electric current is 2el c and the capacitive (or "charging" ) energy is Eq = (2e) 2 /2C, 
where C is the electrostatic capacitance. For trapped BEC's, Ec can be obtained to a good accuracy from the 
Thomas- Fermi (TF) calculation (9) of the chemical potential [see eq. (23)] 



When both Ec and ksT are Ej, the Josephson pendulum Hamiltonian can be approximated as a harmonic 
oscillator, 

H{y,N)~ ±E Jlp 2 + ±EcN 2 , (4) 

whose frequency 



^EjEc/h (5) 



is called the Josephson plasma frequency. 

Given the intimate relationship between eqs. (1-2) and eq. (3), a derivation of the pendulum Hamiltonian (in 
particular of the — cos tp dependence of the energy on the relative phase) may be regarded as equivalent to a derivation 
of the Josephson effect. At low temperatures, the collective dynamics of an inhomogeneous BEC is well described by 
the Gross-Pitaevskii (GP) Hamiltonian Gfl] 
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= I dr(^|-|VvI/(r)| 2 + y 0Xt (r)|vI/(r)| 2 + ||*(r)| 4 ), (6) 

where g — 4:7rh 2 a/m (a is the s-wave scattering length). Neglecting depletion H] , the normalization can be taken as 
j drp = Nt, being Nt the total number of atoms. The extremum condition 8(H — poNt)/S^* = (where /io is the 
equilibrium chemical potential) leads to the well known GP equation flQ|| . 



II. PHASE RIGIDITY AND BROKEN GAUGE SYMMETRY 



The assumption of a uniform phase within each Bose condensate is implicit in most discussions of the Josephson 
effect. It seems therefore convenient to check explicitly that this is a physically meaningful approximation. To prove 
it, let us estimate the energy of a one-radian fluctuation of the phase across the condensate near equilibrium in a 
single spherical harmonic well [see eq. @], 

where R — ao(15oA r r/ao) 1 ^ 5 is the cloud radius estimated within the TF approximation |1, loq is the harmonic 
oscillator frequency of the well, and ao = (S/muo) 1 ' 2 is the oscillator length. For the last approximate equality we 
have used typical parameters a = 5 nm and ao = cm. The characteristic temperature of such an oscillation can 
be as big as 10 /iK 0. Since temperatures as low as 100 nK can be reached in current BEC setups, it is acceptable 
to view the phase of a condensate within one well as internally frozen and rigid. It is clear from (^) that spatial 
fluctuations of the phase are energetically more costly the higher the value of the condensate density p. For this same 
reason, in the interestitial region between two wells, where p(r) decreases appreciably, spatial phase variations require 
less energy and are thus easier to activate at low temperatures. Here lies the essence of the low energy Josephson 
dynamics. At low temperatures, the picture emerges of two wells, each with an internally uniform phase, whose 
relative phase 

<p = J dl-V0, (8) 

induced by local variations of 4>(r) in the interstitial region, is the only active dynamical variable. 

From the theory of phase transitions, we know that some type of "rigidity" is generally exhibited in the ordered 
phase below the critical temperature JT2"| |. In the case of superconductors, superfluids, and Bose- Einstein condensates 
(BEC's), one may speak of "phase rigidity" because, given an electromagnetic gauge (or its equivalent for superfluids 
and BEC's), the phase is determined everywhere once it has been fixed at a given point. Alternatively, since gauge 
and phase are intimately connected in quantum mechanics, we can say that, given a choice of phase, we have lost the 
freedom to choose the gauge, hence the term "gauge symmetry breaking" . 

Long range phase coherence can be destroyed by thermal or quantum fluctuations. A simple picture of how this 
occurs can be developed by analyzing the Josephson Hamiltonian (^|), which may be viewed as a two-site reduction 
of the more general energy functional @. In equilibrium, good phase coherence between sites 1 and 2 is obtained in 
the harmonic limit (Q). Then it is appropriate to speak of global phase coherence or rigidity in the combined double 
BEC system. Thermal fluctuations randomize the phase for fc^T > Ej. To see how quantum fluctuations can destroy 
the two-site phase coherence, we note that, since iV and if are canonically conjugate variables, N can be represented 
as N — id/dip in a quantum description. As a consequence, the capacitive term in (^) plays the role of the kinetic 
energy for the phase variable. At zero temperature, quantum fluctuations can destroy phase rigidity between sites 
1 and 2 if the interactions are sufficiently strong (Ec S> Ej). In this limit, the relative particle number becomes 
a good quantum number [see eq. (JsJ) ] and the two condensates behave as essentially independent. Because of the 
weak connection, it is much easier to destroy the relative phase coherence between condensates 1 and 2, than the 
internal phase rigidity within each condensate. We have already proved that internal coherence cannot be destroyed 
by thermal fluctuations at the experimental temperatures. On the other hand, long range coherence within each 
condensate is preserved from destruction by quantum fluctuations thanks to the three-dimensional character of the 
problem, as stated by the Mermin- Wagner theorem in quantum statistical mechanics. 



III. EXTERNAL JOSEPHSON EFFECT 
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A. Semiclassical description 



Let us assume that V^ x t(r) is such that our system has the structure of two weakly connected condensates 1 and 
2. We assume that the condensates are confined within spherical harmonic wells of the same frequency loq. First we 
wish to analyse the semiclassical dynamics. Then, N and ip can be treated as simultaneously well-defined, and we 
may write for the wave function the ansatz Jl3| 

(r, t) ~ *i(r; N T /2 - N(t)) + e^ w * 2 (r; N T /2 + N(t)), (9) 

where "^(r; n) is the (real) equilibrium wave function for the isolated well i containing n bosons. It is straightforward 
to show that, to lowest order in the overlap integrals J ^1^2, the energy functional for takes the form 

H(tp,N) ~E B (N) + Ej(N){l-cos<p), (10) 

where Eb(N) is the bulk energy of the two isolated wells with TV transferred atoms, and Ej(N) is the Josephson 
coupling energy. 

In the TF limit , the bulk energy Eb (N) is mostly due to interactions, and it may be expanded as 

E B (N) ~ E B (0) + Mo(^f> 2 + ^o"(^f> 4 - (11) 

Since /xo ~ N^/ 5 J^], the ratio between the third and second terms in the expansion is 0.32 N 2 /N^ which means that 
the last term can be neglected in a wide range of situations. To avoid complications stemming from possible resonances 
between Josephson oscillations (see below) and intrawell excitations, we require p,o(Nt/2 + N) — pq(Nt/2 — N) -C Tujq, 
where we use the result that the first normal mode of a spherical well lies approximately at Hluq above the ground 
state [0,^5|. This condition is realized when N/N? <C 4.6A r T 2 ^ 5 for typical parameters. This upper bound is of order 
2 — 10% for Nt ~ 10 4 — 10 6 . 



B. Josephson coupling energy in one dimension 

The expression for the coupling energy Ej(N) in terms of ^ and 9 2 implicit in ( |ic| ) is rather complicated and 
difficult to handle, so it is desirable to find a simpler expression. Let us assume that the (1 — cosp) dependence for 
the energy of a Josephson link has already been proved by e.g. the arguments in the previous subsection, and focus 
on the particular case ip <C 1, when \I>(r) is very close to the ground state wave function \to( r ) with f = 0. In such a 
case, the integrand |V\1/| 2 in (|^) is mostly due to phase variation. In one dimension, we can then write 

E[4>] = f- f dxp(x) \<j>\x)\ 2 ~ l -Ejp\ (12) 
2m J 1 2 

where the phase term in (|j|) has been approximated quadratically, and where the integration extends between two 
points beyond which the phase is practically uniform. 

From (1), we know that a supercurrent must flow between the two condensates whose value is 

/ ~ (Ej/h)<p. (13) 

On the other hand, particle number conservation requires 

I = (h/m)p(x)<j)'(x) = cnst. (14) 

Inserting the resulting expression for <f>'(x) into (|l^), and using (|l3|), one obtains a functional expression for the 
Josephson coupling energy in terms of the ground state superfluid density: 



*2 r 

m 



dx 



Po(x) 



(15) 



where we have used the fact that, for ip -C 1 , p ~ p 
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C. Three dimensions: Electrostatic analogy 



We assume again that the ( 1 — cos ip) dependence has already been proved and consider the case Analogously 
to (|l2|), we can write 



2m 



drp(r)|V0(r)| ; 



(16) 



where the integration extends over the region exterior to the condensates (the results are quite independent on the 
precise location of the condensate borders). The phase within the condensates is assumed to be uniform. The only 
way for eq. ( |l6| ) to have such a dependence on the total phase difference ip and not on the details of V</J is that the 
condition 5E[(p]/6<p(r) — (fi. 2 /2m)V • (pV0) = is satisfied, something which is guaranteed by current conservation 
[(h/m)pV(f> is the particle current density]. The alert reader may have noted that, except for trivial factors, this is 
the electrostatic equation for the electric displacement vector D = p\7<p in a medium with a nonuniform dielectric 
constant p(r). Boundary conditions for <f>(r) are given by its value at the borders of each condensate (ipi and <p%), 
which act as conductors in this analogy. We have a system of two conductors held at a potential difference tp and a 
dielectric medium surrounding them. Then E[<f>] is essentially the energy of this capacitor. Potential theory tells us 
that E[<p] ~ tp 2 , with Ej in ( |l6|) playing the role of the mutual capacitance of the two conductors in the presence of 
such a dielectric. 

Standard variational arguments can be invoked to prove that, for parallel plate boundary conditions [n3[, 



dxdy 



/i dzp{x, y, z)~ 



mE j 
< — J- < 



dz 



J dxdyp(x,y, z) 



(17) 



(here z is the longitudinal coordinate connecting the centers of the two wells, and x, y are the transverse coordinates). 
The lower bound is obtained by removing the positive term (d<j) / dx) 2 + (dtp / dy) 2 from the energy functional (|l6|), while 
the upper bound is derived by taking ip independent of x, y. If the equilibrium density can be factorized in the region 
controlling the capacitance (a not very restrictive condition in practice; see above), we can write po( r ) — /(^j y)ff( z ) 
and the two bounds become identical. Taking r = at the middle point of the double well configuration, and choosing 
x = y = to be the line connecting the two well centers, we can finally write pM 



-2 r 



Ej 



A- 



dz 



Po(0,0,z) 



(18) 



being A — /(0, 0) -1 / dxdyf(x,y) an effective area. The factorization condition is much more general that it might 
appear at first sight, since it only has to be satisfied in the interstitial region near the top of the barrier giving 
the dominant contribution to Ej. There T4 x t(r) has a saddle-point and thus it can be decoupled in longitudinal 
and transverse coordinates. Since interactions are negligible in the bottle neck region connecting the two wells, we 
conclude that the wave function must be factorizable. 



D. WKB calculation 



In order to proceed further, we need an estimate of po = \^o\ 2 in the region of interest. Here we reproduce sketchily 
the WKB calculation of ref. |i3[| . For two identical wells in equilibrium, the ground state wave function is symmetric 
in z. Focussing on the x = y = line giving the dominant contribution, we write 



*o(«) = 



B 



: cosh 



dz'p(z') 



(19) 



where p(z ) = [2m(T4 xt (0, 0, z) — ^o)] 1 ^ 2 and B is a constant to be determined from the work of ref. [^6|. Introducing 
(H) into © we obtain 



Ej 



hA\B\ 



2 tanh ( - 



(20) 



where S = J 2 p(z)dz/h is the dimensionless instanton action. One finally gets |l5| ] 
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e S (N T \ 1/3 /l5a\- 2/3 ^ 



Ej ~ tanh(S/2) ' ~ (21) 



The iY,J/ 3 dependence can be understood qualitatively by noting ]l3||l6|] that A ~ i? 2 / 3 and B ~ i? 1 / 2 . Since, in 



the TF approximation R ~ N^/ 5 , one obtains -Bj ~ N^/ 3 . Knowing that the critical temperature satisfies 
1cbT c ~ N^/ 3 hu>Q JlTfl , we may rewrite ( ^i"| ) as (for large 5 and typical parameters) 

£j ~ k B T c e~ s (22) 

which is reminiscent of the Ambegaokar-Baratoff formula for superconductors Q if one notes that e~ s is the trans- 
mission amplitude for a particle to traverse the barrier. 



E. Josephson plasma frequency 

From (3) and flll|), we have Ec = 2/i' (Nt/2), which, in the TF approximation || reads 

4 /7V T \~ 3/5 / 15a \ 2/5 Tiuj 



From (|), (|T|) and (§|), we obtain 



U -t- (23) 



( 2a \ 2/15 e- S/2 
W " * (,15^ J Vtanh(g/2) W °' (24) 



Since, usually, loq/2-k ~ 10 — 100 Hz, we conclude that iu,jp/2ir < 10 Hz. 

The ratio Ej /Ec is a good measure of the classical character of the relative phase tp. From ( f2l|) and (p3|), we find 

£j / 2a \ 1/15 iVrao e~ s 

E c ~\l5aN T ) a tanh(5/2) ' 1 J 

By varying S, the system can be driven from the classical regime (Ej 3> Ec) to the strong quantum limit (Ej <C -Ec)- 
However, in the latter case, quantum effects are only important if we operate at ultralow temperatures fc^T < .Ec 1 . 
From (|2q), one concludes that, for typical numbers, S must be roughly > 10 for quantum fluctuations to be important. 



IV. EFFECT OF DISSIPATION 



So far, we have described the dynamics of a conservative system. In real life, however, one should expect a certain 
amount of damping. The most obvious source of such damping is the incoherent exchange of normal atoms, and a 
quantitative discussion requires a generalization of our results to nonzero temperature. It is possible to perform a 
qualitative exploration of the effect of damping based on a simple transport model []l3|| . We can expect thermally ex- 
cited normal atoms to follow the usual laws of dissipative currents. Specifically, they must give an Ohmic contribution 
to the current, 

/„ = -Gn, (26) 

where /i is a spontaneous or induced chemical potential difference (for simplicity, we assume equilibrium between 
condensate and thermal cloud within each well). As a result, the first Josephson equation is modified to read 

d Ej 

—N= -±sm<p + I n , 27 
at n 

or, for small fluctuations [using (2) and (pc|)], 

ip + GE C (f> + uj 2 j P p = 0. (28) 
The classical harmonic oscillator described by eq. ( p8| ) is underdamped if 
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It only rests to calculate G for a specific system. A simple estimate can be made for the case in which the chemical 
potential lies near the top of the barrier, at a distance V± ~ fouio/2, so that S ~ 1 — 5. In this regime, we can still 
expect the WKB formula for Ej derived above to yield a reasonable approximation. However, if k B T ^> huj /2 ~ V\, 
the thermal cloud lies mostly above the barrier and a radically different approach is needed to study its transport 
properties. For simplicity, we introduce the drastic approximation that particles impinging on the barrier with energy 
E are transmitted with probability one if E > Vq, and zero if E < Vq. Then, the flow of normal atoms due to a 
fluctuation in [i is only limited by the "contact resistance" , a concept taken from ballistic transport in nanostructures 
jL8[. We may write G = N c h/h, where N c h is an effective number of available transmissive channels. Within a 
continuum approximation, and assuming that only transverse channels with a minimum energy between /xo and 
Mo + k B T are populated, we find (taking /i <C k B T) N c h = mA n k B T/2irh 2 , where A n is a mean transverse contact 
area seen by normal bosonic particles p3[ . Approximating p"H ] A n ~ 2irk B T/mujQ, we find ~ (fc^T /Tiuiq) 2 , and 
hence 

a ~2irNZ /15 (huj Q /k B T) 2 e- s / 2 . (30) 

Interestingly, the same result is formally obtained in the high barrier limit, where normal atom tunneling also lies in 
the WKB regime |0|. This gives additional confidence in the adequacy of the estimate (|30"|). Taking k B T = lOftwo, 
we find a ~ 0.38 for N ~ 10 4 , and 3.3 for N ~ 10 6 , if S ~ 5. The equivalent numbers for S ~ 1 are 2.8 and 24. The 
conclusion is that coherent Josephson dynamics can be observed in current atomic Bose condensates if the barrier is 
low. Underdamped dynamics can be further favored by decreasing T and increasing N . 



V. INTERNAL JOSEPHSON EFFECT 



It is possible to couple different hyperfine states of an atom (say \Fm > and \F'm! >) through laser light. In 
particular, by applying a laser pulse, one may force an atom initially in state \Fm > to evolve in a controllable way 
towards another state \F'm' >. If such a pulse is applied to a macroscopic condensate, then it is the whole ensemble 
of atoms what evolves coherently. Thus it is possible to prepare a condensate in which each atom is in the same 
coherent superposition of the two states. Such an experiment has been realized by Hall et al. Jl9t with the |1, — 1 > 
and |2, 1 > states of 87 Rb (hereafter, \A > and \B >). A natural variant of this experiment consists in using laser 
light which also induces the opposite rotation {e.g. linearly polarized light), in such a way that the two states are 
coherently connected in both directions, very much like Cooper pairs in a Josephson junction can tunnel from left to 
right and vice versa. In such a setup, it will be possible to study the novel and, potentially, extremely rich "internal" 
Josephson effect, which bears some analogies with the physics of superfluid 3 He-A, as noted by Leggett [pp| . 

The analysis of the spatial (external) Josephson effect given in section III is based on the knowledge of I4 x t(r), 
where r we can be varied continuosly between wells. A similar study of the internal Josephson effect would not be 
practical. Fortunately, all is really needed is the value of the matrix element connecting states \A > and \B >, which 
can be known by other means. Then, it is most convenient to employ a two-site description of the Josephson link. 
We write HQ 

H = -^(a^b + &t ) + |[( a t a ) 2 + (6+6) 2 ], (31) 

where a' (&') creates one atom in state \A >(\B >), u accounts for the interactions, and ojr is the Rabi frequency. For 
simplicity, we assume that the two atomic energies, as well as the intraspecies interactions, are equal, and neglect the 
interspecies interaction. Particle number conservation requires a^a + b^b = Nt = 2Ao. 

To establish the connection with the pendulum Hamiltonian, we note that, since particle eigenstates admit a 
representation $jv(v) = (27r)~ 1//2 exp(— iNip), we may write for (3) 

<N+l\H\N>=-Ej/2. (32) 

A similar analysis for the two-site Hamiltonian (^Tj) yields 

< N + 1\H\N > = -(Rwh/2) <N A + 1,N B - l\a%\N A , N B > 

= -{hu R /2)^{N A + l)N B . (33) 



G 



Noting that Na,b — Nq ± N , and identifying ([32]) and (|33|), we arrive at 

Ej(N) = huj R [N (N + 1) - N(N + 1)] 1/2 . (34) 
In the limit 1 < N < jV 0) eq. @ becomes 

Ej = NoHwr, (35) 

which is a clear manifestation of the phenomenon of Bosonic amplification. The identification of (|3l]) with the 
pendulum Hamiltonian is completed by noting that (apart from constant energy shifts) the interaction term can be 
written uN 2 , i.e. u = Ec/2- 



VI. CROSSOVER FROM JOSEPHSON TO RABI DYNAMICS 



A clear limitation of the standard pendulum Hamiltonian (3) is that, in the non-interacting limit Ec — * 0, the 
dynamics is suppressed; in particular, the Josephson plasma frequency (|^) for small collective oscillations vanishes. 
This contradicts our physical notion that noninteracting atoms should indeed exhibit some dynamics. For instance, 
under the effect of Hamiltonian (|3l|), an atom initially prepared in state \A > will undergo Rabi oscillations between 
states > and \B > with frequency uir. This effect is clearly beyond the scope of a rigid (with TV-independent Ej) 
pendulum model such as (3). However, it should be describable by (j3ll), which has a well-defined non-interacting limit. 
Thus, it seems of interest to develop a unified description which views collective Josephson behavior and single atom 
Rabi oscillations as particular cases of a more general dynamics. To that end, let us rewrite (|34|) in a semiclassical 
form, 



Ej(N) = hu; R ^N$-N 2 . (36) 
The classical dynamics of the resulting non-rigid pendulum Hamiltonian, 



H = -N huj R S Jl- N 2 /N 2 cos cp+{E c /2)N 2 , (37) 



has been analyzed by Smerzi et al. |23f| . In the limit of small <p and N, (|37]) acquires the harmonic form 
with a natural oscillation frequency 



H = N hcj R ^+ [Ec + j^ J — , (38) 



uj 2 



J .JP ' ^R 

.2 



(N Ec/h)toR + to R. (39) 



It is thus clear that, for a given interaction strength, the double BEC system can be driven continuously from 
the Josephson to the Rabi regime by varying uift, something feasible with current laser technology. Noting that 
Ec ~ L^o/No, the Josephson limit corresponds to (Iq 3> lur. 

The difference between Josephson and Rabi dynamics is more dramatic when one considers large values of N, 
comparable to A^o- When Ec = 0, it is easy to see that (j3?j) always yields oscillations of N(t) around zero with 
frequency u>r. By contrast, when Ec ^ 0, the frequency depends on the initial value N(0) and, when the ratio 
A = NqEc /hujR exceeds a critical value A c , the oscillations average to a nonzero number |£Ug3[ . This number A c 
depends on Af(0) and ^(0). In particular, A c = 2 when A^(0) takes its maximum value A^o, and A c — > oo when 
N(0) — - > 0, i.e. in the small amplitude limit, there are always oscillations of zero average and frequency given by 
(p3). The phenomenon by which, due to macroscopic quantum coherence, it is possible (for A > A c ) to create a large, 
long lived particle number imbalance between two connected BEC's, is called quantum self-trapping and it is 

equivalent to the nonzero average angular momentum of a mechanical pendulum undergoing complete turns. Like the 
pendulum motion, the average particle imbalance can decay because of dissipation. Quantum self-trapping cannot be 
realized in superconducting Josephson junctions because, the chemical potential differences would have to be greater 
than the superconducting gap |p3[ , thus allowing quasiparticles to intervene and complicate the dynamics. This is not 
an important limitation in the case of BEC's. There exists a close and well-studied analog of quantum self-trapping 
in the longitudinal NMR of 3 He- A pf. 
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Eq. (|37]) describes in a unified way Josephson and Rabi (Ec — 0) dynamics. This is analogous to how eq. @ can 
yield both the Gross-Pitaevskii and the Schrodinger (g — 0) equations for the wave function of a many boson system. 
The crossover between collective Josephson and individual Rabi dynamics cannot be studied in superconductors 
and supcrfluids, because there interactions are never completely negligible. It is a nice feature of Bose-Einstein 
condensation that it will allow us to study the crossover between these two qualitatively different dynamical regimes 
in an elegant fashion. 

VII. PHASE DYNAMICS OF SEPARATE CONDENSATES 

The class of Josephson links describable by the harmonic approximation (||) includes many of the foreseeable 
connections that will be realized between BEC's. At equilibrium, the phase for such systems is narrowly peaked 
around ip — 0. Being canonically conjugate variables, N and ip satisfy the uncertainty relations ANAip > 1/2. In 
general, their m.r.s. values are 

Ap=(E l /E J ) 1/2 AN=(E l /E c ) 1/ \ (40) 

being Ei = (Tiuijp/2) coth.(hujjp /2ksT) the average energy of the oscillator. If Ec and fc^T <C Ej, we have Aip <C 1. 
It is then of interest to analyze the phase dynamics of two BEC's which, having formed an equilibrium Josephson 
link, are suddenly disconnected [Ej is made zero). Such a study may allow us to understand better the robustness 
of the phase definition when the support of the Josephson coupling is no longer present. One may naively think that 
two separate BEC's have a well-defined relative particle number and that the phase is necessarily random. However, 
if the starting point is that of a well-defined phase, the situation is much less clear, since phase randomization needs a 
finite time that may be sensitive to the details of the environment. After all, we are all familiar with physical variables 
(such as the position of macroscopic objects) which are permanently well-defined and yet are not constants of motion. 
There are at least three physical factors which favor this long term definition of non-conserved quantities: (i) large 
inertia (small effective Ec), (ii) high potential barriers in the effective potential landscape, and (iii) damping due to 
the environment. For the phase of a Josephson junction, at least the first two ingredients (and some times the three 
of them; see section IV) are present. Here we wish to understand what is the fate of the phase when the second factor 
(the Josephson coupling) is removed. 

The randomization of the phase after suppression of the Josephson coupling was investigated in refs. Q and [g5| for 
the case of superconductors and superfluids. The possibility of performing interference experiments between BEC's 
p6| , [l9f has revived the interest in the study of this |^7) and related problems . 

A. Ballistic randomization 

In the absence of external perturbations, the evolution of the wave function at t > is that of an initially narrow 
wave packet on a circle. Assuming for the moment that if is an extended variable, the phase uncertainty at long times 
is 

Atp(t) ~ y/EiEct/h = ANE c t/h (41) 
The ballistic time tb is defined as the time it takes Aip to be of order 2ir. We get 

T B = h/y/EiEc. (42) 

A typical number for superconductors can be tb ~ 10~ 8 s, while, for superfluids, the equivalent time may be of the 
order of minutes or hours, depending on whether the capacitive forces are due to compressibility (fixed volume) or 
gravity p^ , p5[ . The randomization of the relative phase between BEC's lies in an intermediate time scale. Using (pS]), 
we find t b ~ 20 - 100 ms (with N ~ 10 4 - 10 6 ) for T = 100 nK. 

B. Recurrent dynamics 

The phase is actually an angular variable for which and 27r are to be identified and this implies that N has integer 
eigenvalues. When Ej — 0, the wave function evolves as 
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$(¥>, t) =< <p\Q(t) >= a N<S>N&) eM-*E c N 2 t/2h), 

N 



(43) 



where a^r =< iV|$(0) >. From p3|), we notice the existence of a recurrence time 

r fi = (44) 

after which the wave packet recombines, $(y>, t^) = $((^,0). The property JSj <C -Ej implies tr ^> tb, so that the 
two time scales are well separated. This observation renders meaning to the apparently oversimplified mechanism of 
ballistic spreading, provided one considers times much smaller than tr. 

A detailed analysis reveals the following picture |^9| : An initially narrow wave packet begins to broaden according 
to (f4l|). After a time of order tb, the wave packet is almost uniformly spread around the circle, apart from fast 
oscillations. At a much later time tr, the wave packet recombines and broadens again in a time scale tb- A typical 
value may be tr ~ 10~ 5 s for a superconductor, and an astronomically large one for a superfluid. Again, BEC's lie 
somewhere in between: We find tr ~ 30 — 500 s for N ~ 10 4 — 10 6 . It must be noted that recurrence would only be 
observed in quite ideal conditions (very low noise, good time resolution, etc.), perhaps not easily realizable. 



C. Diffusive randomization 



So far we have neglected the possible effect of dissipation due to the environment. In refs. |25|]29| ], a source of noise 
was studied which, when really present, can be expected to dominate the dynamics, namely, the exchange of particles 
with a normal reservoir. In BEC's, such a reservoir could be provided by the thermal cloud. 

Following ref. p5[ , we base our study on a useful analogy which can save us a good amount of work. If the contact 
between the condensate and the reservoir is Ohmic, i.e. the rate of particle exchange behaves as 

SN = -TSfi. (45) 

where Sfj, is a possible chemical potential difference between condensate and thermal cloud, then our problem is 
formally identical to that of a Brownian particle (BP). The phase plays the role of the position and the capacitance is 
analogous to the mass. The particle number N is equivalent to the momentum p (both are conserved quantities that 
are exchanged with the bath), and the chemical potential difference \i ~ <fi is the the analog of the particle velocity v 
(the relative position of two BP's behaves like another BP). The friction coefficient rj, defined as p — —rjv (where v is 
a deviation of the velocity from its zero equilibrium value), is equivalent to the coefficient T. 

Once the analogy with the BP has been established, we may borrow from the classical theory |3(J and write for the 
evolution of the phase uncertainty 



w h 2 r 



i-i(l-c-rt) 
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(46) 



where 7 = TEc- At short times (7* <C 1), we reproduce the ballistic randomization given by pi]). By contrast, at 
long times (jt 3> 1), the randomization becomes diffusive, 

Atp(t) ~ (2k B Tt/h 2 T) 1/2 . (47) 

Eq. ( ]47| ) is a classical result. However, thanks to the work of ref. |H| on a quantum BP with Ohmic dissipation, we 
know that the diffusive law ( [47|) applies at long times t 3> Ti/ksT ^S> 7 _1 , provided A</?(0) <C 1 
Analogously to the ballistic case, it is possible to define a time for diffusive randomization 

t d = h 2 T/2k B T (48) 

For superconductors, td may be of the order of minutes, a surprisingly long time if we compare it with the estimates 
for tr. This is a very important effect. The large disparity between the two time scales may be understood by 
noting that, while the fast ballistic randomization is due to the large value of Ec (small electrostatic capacitance), the 
diffusive dynamics (which takes over at a time 7 -1 <C tb) is essentially independent of it and is determined instead 
by the coupling to the normal reservoir. We arrive at the apparently paradoxical result that, upon suppression of the 
Josephson contact, dissipation may help to preserve phase coherence. 







D. Common vs. independent thermal clouds 

It is important to emphasize that the diffusive mechanism described above only applies if the normal particle 
reservoir is common to both condensates. In VII. C, we have considered a situation in which, although the two 
condensates have been separated, the two thermal clouds remain connected to the point of behaving as a single one. 
In superconductors, an equivalent arrangement is not fundamentally difficult (the reservoir would be a normal metal). 
In the case of BEC's, we may have a common thermal cloud when the increase in the barrier height is large enough 
to render the two condensates unconnected, but still sufficiently small to allow for an essentially single thermal cloud. 

If the two reservoirs are independent, then, following the BP analogy (in which the finite masses of each of the 
separate baths would cause a ballistic spread of the uncertainty in the relative position), the finite effective value of 
the total Ec for each one of the condensate plus reservoir system will give rise to a ballistic randomization of the 
relative phase. This seems to be the case in interference experiments between different hyperfine states [ fl9| . 

The crucial difference between ballistic and diffusive spreading of the phase is that, in the former case, there is no 
restoring force for fluctuations in N . Therefore, any initial deviation AN gives rise to a permanent A/z which causes 
a ballistic spread of the phase [see (ftl|)1. On the contrary, if randomization is diffusive, an initial departure from 
equilibrium is quickly compensated and, in the average, N(t) and ji(t) change sign many times, in such a way that 
Aip(t) grows slowly [see J47|)]. 

VIII. CONCLUSIONS 

I have presented an overview of the physics of the Josephson effect betwee Bose condensed systems, with emphasis on 
the recently achieved BEC's in trapped alkali gases. I have mostly focussed on those physical phenomena that are likely 
to be observed only (or more easily) in these novel systems. So I have omitted the discussion of problems (such as e.g. 
steady particle flow under the action of an alternating chemical potential) which may be viewed as straightforward 
applications of well known Josephson physics I have tried to underline the potential richness of the physics 
displayed by weakly connected BEC's. One may have an external (spatial) and an internal (hyperfine) Josephson 
effect. It seems possible to explore the crossover between collective Josephson behavior and independent boson 
Rabi dynamics. Finally, the observation of fascinating phenomena such as quantum self-trapping and macroscopic 
interference between separate Bose condensates seems also within reach of the emerging BEC technology. Everything 
indicates that the experimental and theoretical study of the Josephson effect between Bose- Einstein condensates will 
lead us to the exploration of most exciting new physics. 
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